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We measure the ferroelectric polarization of BiFeOs films down to 3.6 nm using low energy electron 
and photoelectron emission microscopy. The measured polarization decays strongly below a critical 
thickness of 5-7 nm predicted by continuous medium theory whereas the tetragonal distortion does 
not change. We resolve this apparent contradiction using first-principles-based effective Hamiltonian 
calculations. In ultrathin films the energetics of near open-circuit electrical boundary conditions, i.e. 
unscreened depolarizing field, drive the system through a phase transition from single out-of-plane 
polarization to nanoscale stripe domains. It gives rise to an average polarization close to zero as 
measured by the electron microscopy whilst maintaining the relatively large tetragonal distortion 
imposed by the nonzero polarization state of each individual domain. 



A major issue for prospective nanoscale, strain- 
engineered fD ferroelectric devices [2| is the decrease of 
the polarization of ultrathin films. The depolariz- 
ing field arising from uncompensated surface charges re- 
duces or even suppresses ferroelectricity below a criti- 
cal thickness [3l [4] . Ferroelectric capacitors for example 
may exhibit a critical thickness E] • Lichtensteiger et 
al [T have shown that the decrease in Pr in PbTiOs 
(PTO) thin films between 20 and 2.4 nm on Nb-doped 
SrTiOs (STO) substrates is concomitant with that of 
the tetragonality (ratio of the out-of-plane to in-plane 
lattice parameter c/a)). On Lao.erSro.asMnOa (LSMO), 
PTO polydomains were formed below 10 nm with high 
tetragonality [8 . The formation of a polydomain state 
has been suggested for SrRu03/Pb(Zr,Ti)03/SrRu03 
capacitors with Pb(Zr,Ti)03 thicknesses below 15 nm [9 . 
Pertsev and Kohlstedt showed the importance of mis- 
fit strain for the critical thickness of the monodomain- 
polydomain stability for PTO and BaTi03 [lO^. Using 
piezoresponse force microscopy (PFM), BiFe03 (BFO) 
films have been shown to remain ferroelectric down to a 
few unit cells pJT[T3] with both the polarization and the 
slope of the piezoresponse hysteresis loop scaling with 
tetragonality. However, PFM is very local and can only 
provide indirect, semiquantitative estimates of the polar- 
ization. Imperfect tip surface contact can contribute to 
polarization suppression via the depolarizing field. Direct 
electrical measurements of the polarization- field (P(E)) 
loop in ultrathin ferroelectric films are a challenge be- 



cause of leakage current for thicknesses below a few tens 
of nm [TTl |T4]. They become impossible in the tun- 
neling regime for ultra-thin films (5 nm or less) which, 
furthermore, is of the same order as the critical thick- 
ness, heff, estimated from Landau- Ginzburg-Devonshire 
(LGD) elastic theory for polarization stability [(131 [15]. 
BFO can accommodate in-plane compressive strain via 
out-of-plane extension and through oxygen octahedron 
rotation about (111) [16], a degree of freedom not avail- 
able in P4mm PTO films. This interplay between strain, 
tetragonality and octahedra rotations leads to an unex- 
pected decrease of Tc with strain, at odds with the vari- 
ation of c/a ratio. Thus the relationship between struc- 
tural parameters and the remnant out-of-plane polariza- 
tion in very thin films remains an open question. 

In this Letter, we have studied the polarization of BFO 
films from 70 to 3.6 nm thick using a combination of x-ray 
diffraction (XRD), mirror electron microscopy (MEM) 
and photoelectron emission microscopy (PEEM). The 
electron microscopy techniques provide full-field imaging 
of the electrostatic potential above the surface and the 
work function whereas the tetragonality is measured by 
XRD. The results are interpreted in the light of a three- 
dimensional (3D) generalization of previously developed 
dead layer model for thin films within the framework of 
continuous medium theory that predicts a fast decrease 
of the polarization when decreasing the thickness. Inter- 
estingly, the extremely low polarization below heff does 
not scale with the tetragonality and is explained using 
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first principles-based effective Hamiltonian calculations 
which show that as a function of screening the films un- 
dergo a phase transition from single to nanoscale stripe 
domains with an overall out-of-plane polarization close 
to zero. 

Bilayers of BFO/LSMO were epitaxially grown on (001)- 
oriented STO substrates by pulsed laser deposition us- 
ing a frequency tripled {hu = 355 nm) Nd-doped yt- 
trium aluminium garnet (Nd:YAG) laser at a frequency 
of 2.5 Hz [U. The 20 nm thick LSMO layer serves as a 
metallic bottom electrode for ferroelectric characteriza- 
tion. XRD measurements on 70 nm to 3.6 nm-thick thin 
films were performed to track the out-of-plane parameter 



and c/a ratio (Fig. la). The c/a increases slightly from 
1.050 for the 70 nm film to 1.053 for 7 nm, then remains 
constant down to 3.6 nm. This contrasts dramatically 
with the behavior of PTO reported in [7 where c/a de- 
creases with thickness. The chemistry of the films was 
measured by X-ray Photoelectron Spectroscopy (XPS). 
Figure [lb] shows spectra from Bi 4f core- level for thickest 
(70 nm) and thinnest (3.6 nm) films. The spectra are 
virtually identical for both films (for intermediate thick- 
nesses, see ^7 ) showing that the chemical state and sto- 
ichiometry do not change. Bi 4f spectra have a thickness 
independent component shifted by 0.6 eV to higher bind- 
ing energy, suggesting that our strained thin films do not 
exhibit the several nanometer thick skin observed on sin- 
gle crystals [18 . C Is spectra show that contamination 
of the BFO surface is similar for every thickness sug- 
gesting a similar contribution to extrinsic screening in all 
films [17]. 

For the 70 nm BFO film, the ferroelectric properties 
were investigated by standard polarization versus elec- 
tric field P(E) loops (Fig. Ic). The piezo-response hys- 



teresis loops are shown in Fig. |ld| They are position 
independent and exhibit similar coercive values as non- 
local P(E) loops, attesting to sample homogeneity. In 
a (001) BFO film P+ and P~ states are the projec- 
tions of (111) polarization along [001]. Poling of micron 
sized domains was performed by applying a d.c. volt- 
age higher than the coercive voltage (inferred from the 
piezoresponse loops) on the tip while the bottom elec- 
trode was grounded. PFM imaging was carried out at 
an excitation frequency of 4-7 kHz and an a.c. voltage 
of 1 V. No morphology change occurred during poling 
as checked by Atomic Force Microscopy. A low-energy 
electron microscope (LEEM, Elmitec GmbH) was used 
to measure the electron kinetic energy of the MEM (re- 
flected electrons)-LEEM (backscattered electrons) tran- 
sition with a spatial resolution of 30 nm. The transition 
energy is a measure of electrostatic potential just above 
the sample surface [19] and depends on polarization and 
the screening of polarization- induced surface-charge [20^ . 
It allows a noncontact estimation of the out-of-plane po- 
larization for tunneling films, otherwise inaccessible to 
standard electrical methods. All experiments were done 
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FIG. 1. (a) Evolution of pseudocubic lattice parameters a 
and c and c/a ratio with thickness, (b) Bi 4f spectra for 
70 and 3.6 nm films showing two components (surface in 
light grey, bulk in dark grey) for each spin-orbit core-level, 
(c) Polarization-voltage and current-voltage hysteresis loop 
of BFO(70 nm)/LSMO(20 nm)//STO(001) (d) Piezoresponse 
hysteresis loop (local measurement under the PFM tip). 




FIG. 2. (a) PFM phase image of written P^/P domains. 
Each domain is 5x20 /xm^, (b) LEEM image for Einc = 
1.40 eV, (c) PEEM image at E - Ep = 4.35 eV. 



at least two days after domain writing to ensure that the 
observed contrast is not due to residual injected charges. 

Figure |2b] shows a typical LEEM image with a field of 
view (FoV) of 33 /im for incident electron energy (Einc) of 
1.40 eV. The observed contrast reproduces well the PFM 
image of Fig. |2a| A full image series across the MEM- 
LEEM transition (E) was acquired by varying Einc from 
-2.0 to 3.0 eV. Figure [3a| displays the electron reflectiv- 
ity curves showing the MEM (high reflectivity) to LEEM 
(low reflectivity) transition for the P+ (brown upwards 
triangles, E = 0.75 eV) and P~ (green downwards trian- 
gles, E = 1.20 eV) domains. Using complementary er- 
ror function fits we obtain MEM-LEEM transition maps 
showing clear contrast in the electrostatic potential above 
the surface between the P+, P~ and unwritten regions 
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FIG. 3. (a) Reflectivity spectra extracted from tfie and 
P~ domains, (b) MEM-LEEM transition map obtained from 
the image series (70 nm thin film). 




FIG. 4. (a) Threshold spectra extracted from and P~ 
domains (b) Work function map obtained from the threshold 
image series (70 nm thin film). 
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FIG. 5. Thickness dependence of (a) A$i? (red squares) and 
AE (black circles), (b) Pz/Pmax calculated from PEEM (red 
squares) and MEM-LEEM (black circles). Red curve is fit to 
PEEM/LEEM data with hefr = 5.6 nm. Blue diamonds are 
Pz/Pmax values used for numerical simulations, (c) screening 
coefficient /3 calculated from experimental Pz/Pmax values. 



(Fig.|3bl). 

The energy filtered PEEM experiments used a Na- 
noESCA X-PEEM (Omicron Nanotechnology GmbH). 
PEEM of the photoemission threshold gives a direct, ac- 
curate (± 20 meV) and nondestructive map of the work 
function [21^ that may depend, for example, on domain 
polarization [22]. Image series were acquired over the 
photoemission threshold region with a mercury lamp ex- 
citation {hu = 4.9 eV). The lateral resolution was es- 
timated to be 200 nm and energy resolution 200 meV. 



Figure [2c] shows a typical PEEM image of the prepoled 
P+ and P~ regions for the 70 nm BFO film. The energy 
contrast between oppositely polarized domains fits the 
PFM image except at the domain boundary where the 
lateral electric field induced by a P+/P~ domain wall 
deflects electrons [23]. We extract the threshold from 
the pixel-by-pixel spectra using a complementary error 
function to model the rising edge of the photoemission 
(Fig. 4a). Figure 4b maps the work function in the P+, 



P and as-grown regions. 

The difference in the MEM-LEEM transition of the 
P+ and P~ regions, AE, varies from 450 meV for the 
70 nm film to 25 meV for the 3.6 nm film and is plotted 
in Fig. 5a (black circles, right axis). The mean work 
function difference measured in PEEM between P+ and 
P" domains, A^^t^ = ^f{P^) - ^f{P~), is plotted as a 
function of thickness in Fig. [5a| (left axis). While A^f is 
300 meV between 70 nm and 7 nm, between 7 and 5 nm 
it drops to 20 meV. 

The polarization charges at the BFO surface are 



screened over a so-called dead layer leading to an inward 
(P+) or outward (P~) surface dipole. By measuring the 
work function (or surface potential) difference between 
two opposite domains, our method allows a direct mea- 
surement of the polarization-induced dipoles since any 
averaged nonferroelectric contribution is canceled. The 
surface dipole difference, hence the surface potential or 
work function difference, is proportional to the difference 
in polarization charges when going from the P+ to the 
P~ domains: 

A$i.,A^(x - {P^ .d-^ - P- .d-) ^2—Pr.d (1) 

Co ^0 

where P+'~ and d'^~ are the polarization and dead layer 
thickness for the upward, downward domains; Pr is the 
average magnitude of the polarization in the two poled 
domains and d is the average dead layer thickness. For 
the sake of generality, one can take into account electronic 
screening via a high-frequency dielectric permittivity, but 
it would still leave a linear relation between polarization 
and A^F, AE. Pz/Pmax, where Pz is the measured out- 
of-plane polarization and Pmax the value for the 70 nm 
film, is plotted as a function of film thickness in Fig. |5b| 
By comparison with Fig. [Taj the drop of average polariza- 
tion between 7 and 5 nm does not result from a decrease 
in the c/a ratio, contrary to PTO thin films 0. Here the 
c/a ratio increases for thinner films and is constant at 
1.054 below 5 nm. If there were no polarization then it 
should be about 1.03. However, PTO is almost fully re- 
laxed whereas BFO is compressively strained. Secondly, 
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in BFO, the polarization deviates appreciably from the 
[001] direction and is the macroscopic average of four 
(111) type distortions. We have therefore generalized 
the ID dead layer LGD model of Bratkovsky and Lev- 
anyuk [15 to the 3D polarization case. It gives the follow- 
ing relation for thickness dependence of polarization [17] : 



A\ B 



h 



(2) 



where heff is the effective thickness below which the 
macroscopic Pz goes to zero, and A, B are fitting pa- 
rameters. A good fit to the data is obtained with heff = 



5.6 nm (see Fig. 5b, red curve), compared with 2.4 nm 
for PTO. 

To understand why the measured polarization sud- 
denly drops in ultrathin strained (001) BFO films, while 
the axial ratio is still very large, we have conducted first- 
principles-based, effective Hamiltonian calculations [24]- 
[26| that take into account free surfaces [24 . We used the 
lattice parameter of the STO substrate for the pseudo- 
cubic in-plane lattice constant of BFO, leading to a mis- 
fit strain of -1.8%, in agreement with the experimental 
value. The calculation includes the local electric dipoles, 
the strain tensor and tilting of the oxygen octahedra. 
The electrical boundary conditions are governed by a co- 
efficient denoted as /3 described in Ref. 27, Practically, P 
can vary between (ideal open-circuit, maximal depolar- 
izing field) and /3 = 1 (ideal short-circuit, fully screened 
depolarizing field) . To determine P for each of our grown 
films we first extract the Pz/Pmax values from a B-spline 
interpolation of the experimental data (Fig. [5b| blue dia- 
monds) and then vary P in the calculations until the pre- 
dicted Pz/Pmax perfectly agrees with the experimentally 
extracted one. Figure [5c] shows the resulting p values. 
P decreases with thickness, indicating that the observed 
decrease of polarization is related to imperfect screening 
of the depolarizing field. The vanishing of the overall z- 
component of the polarization (which occurs experimen- 
tally for thicknesses lower than 5.6 nm, see Fig. 5b) is 



associated with values of P lower than 0.4 (see Fig. 5c). 

To understand what happens for these P values, 
we performed additional first-principles-based effective 
Hamiltonian calculations on a single 20x20x20 super- 
cell (i.e. with a thickness of 8 nm) allowing P to vary. 
This supercell was chosen because around 8 nm the po- 
larization is very sensitive to the thickness (Fig.|5b|). The 
results are shown in Fig. [6] At a critical value of P of 
0.275 ± 0.025 the BFO supercell goes from a phase with 
a uniform out-of-plane polarization to a stripe domain 
phase with a vanishing overall out-of-plane polarization. 
Fig. [6a] displays the energy of these two phases as a func- 
tion of p. The monodomain phase is energetically more 
favorable than the stripe nanodomains for P above 0.30 
and less for smaller P values. The predicted evolution 
of the c/a ratio, and of the overall Pz/Pmax, with P for 
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FIG. 6. (a) (top panel) Slab energy with P showing a tran- 
sition from single domain (out-of-plane polarization, green 
triangles) to stripe domains (no total out-of-plane polariza- 
tion, red squares) below pcrit = 0.30, (bottom panel) theo- 
retical c/a ratio with P for the two phases, (b) Evolution of 
Pz/Pmax with P, inset: domain morphology (upper left, stripe 
domains; bottom right, single domain). All these data are for 
a 20x20x20 slab, i.e. a film with a thickness of 8 nm. 



single and stripe domain phases are shown in Figs. 6a 
and [6b] respectively. Interestingly, a continuous ferro- 
electric to paraelectric transition would lead to a large 



monotonic decrease of tetragonality (Fig. 6a, green tri- 
angles), which we do not measure below heff. Rather, 
the transformation from ferroelectric monodomains to 
nanostripe domains leads to a (large) c/a similar to the 
one associated with short-circuit-like conditions (i.e. for 
which P is close to 1). Such results are consistent with 
our experimental findings that c/a does not vary between 
70 nm and 3.6 nm, and explains that such insensitivity 
to strain is likely due to the formation of nanostripe do- 
mains. The single to stripe domain transition explains 
the loss of contrast in electronic microscopy observed in 
LEEM and PEEM contrast between 7 and 5 nm, because 
these regions do not possess any overall z-component of 
the polarization. The stripes have a typical dimension of 
a few nanometers, which is below the lateral resolution 
of our experiments (The top left inset of Fig. 6b shows 
the morphology of these domains). However, stripe do- 
mains in BFO thin films close to the heff value have been 
observed by PFM [28 . For such thin films, one might 
also ask to what extent the screening at the LSMO/BFO 
interface affects the measured polarization. Transmission 
electron microscopy of the interface between LSMO and 
a 3.2 nm BFO film suggests that the first three BFO unit 
cells are screened by the interface charge [29]. This also 
fits nicely with our experimental observation of an abrupt 
decrease in polarization starting at 7 nm, 1.4 nm above 
the calculated heff. 

In summary, we have measured the polarization in 
ultrathin strained (001) BFO films using PEEM and 
LEEM. The polarization drops abruptly below a criti- 
cal thickness hcrit whereas the tetragonality has a high 
constant value. A first-principles-based effective Hamil- 
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tonian approach suggests that BFO exhibits a first order 
phase transition to stripe domains at hcrit = 5.6 nm, cor- 
responding to a screening factor, /3, below 0.35. This 
model fits the experimental measurement of the average 
polarization and the c/a ratio. 
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SUPPLEMENTARY MATERIALS 



X-Ray PhotoEmission Spectroscopy for Every Film Thickness 



X-ray PhotoEmission Spectroscopy (XPS) was carried out using a Kratos Ultra DLD with monochromatic Al Ka 
(1486.7 eV). The analyzer pass energy of 20 eV gave an overall energy resolution (photons and spectrometer) of 

0. 35 eV. The sample is at floating potential and a charge compensation system was used. The binding energy scale 
was calibrated using a clean gold surface and the Au 4f7/2 line at 84.0 eV as a reference. A take-off angle of 90°, 

1. e., normal emission, was used for all spectra presented. The XPS spectra show that the chemical environment was 
identical within 1% (see 7a and 7b) for every film. Krug et al. [20] pointed out the importance of adsorbates on LEEM 
and PEEM measurements. Figure 7c shows that surface contamination is similar in 3.6 (low contrast) and 70 nm 
(high contrast) thin films indicating strongly that the disappearance of ferroelectric contrast is not due to differential 
contamination. Moreover, the 5 nm film has the lowest carbonates concentration and still shows weak ferroelectric 
contrast in LEEM/PEEM experiments (see Table 




168 164 160 156 728 718 708 288 286 284 282 

Binding Energy (eV) Binding Energy (eV) Binding Energy (eV) 

(a) Bi 4f (b) Fe 2p (c) C Is 



FIG. 7. XPS spectra of BiFeOs thin films for every thickness. 



Thickness (nm) 


Icis 1 ^Bi4f 


3.6 


2.22 


5.0 


1.49 


7.0 


3.30 


20 


1.40 


70 


2.20 



TABLE L C Is to Bi 4f ratio calculated from XPS spectra 
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PFM/PEEM/LEEM Data for Every Film Thickness 



In addition to the 70 thin film data presented in the main manuscript, Fig. [8] displays the Piezoresponse Force 
Microscopy (PFM) images for every thickness showing the films have been successfully poled. For the thinnest film 
(3.6 nm), Fig. [9] shows the PFM phase loop as a function of d.c. bias. During writing process, we applied an electric 
field at least twice as high as the coercive field. 




(a) 70 nm (b) 20 nm (c) 7 nm (d) 5 nm (e) 3.6 nm 

FIG. 8. PFM images for every thickness. 




-2-10 1 2 
Voltage (V) 



FIG. 9. PFM phase loop for the 3.6 nm thin film. 



Fig. 10 and 11 show work function maps and MEM-LEEM transition maps for every thickness. They show clearly 
the different behavior between the 70, 20 nm thin films (high work function difference between P+ and P- poled 
regions) and the 5, 3.6 nm thin films (near-zero difference in work function). The intermediate case of 7 nm shows 
clear contrast (left part of Fig. 10c for instance). At the same time, clear contrast is visible within the P+ and P- 
(especially in the MEM-LEEM images) poled domains. Although the lateral resolution does not allow visualization of 
the nanodomains found by the numerical simulations, this additional contrast may be indirect evidence of the switching 
process, indicating that the formation of the stripe domains is not simultaneous across the full width of the poled 
domain. Oxygen vacancies can pin domain walls in, for example, PTO ^ and (Bio. gsPro. 15) (Feo.gsMno. 05)03 |3T] . 
The presence of such defects could therefore act as a nucleation center for the initial stripe domain walls. The 
formation of the stripe domains would then proceed outwards from the initial domain walls; indeed there is evidence 
in both PEEM and MEM-LEEM images for a much finer striped structure within each poled domain. A detailed 
study of the stripe formation is beyond the present manuscript. We believe that this illustrates single domain vs stripe 
domains competition around the transition thickness, likely due to slightly different boundary conditions and/or to 
the presence of defects in the film. It will be the subject of future work. 



Threshold Spectra Using X-ray Source 



Photoemission at threshold shows a cut-off energy below which electrons cannot escape the surface. It is often 
assumed that only secondary electrons contribute to the threshold spectra and that the complementary error function 
(erfc) is the correct function to deduce the work function from the rising edge of the photoemission threshold. However, 
the emission spectrum of the Hg light source peaks at only 4.9 eV. With such low photon energy, direct transitions 
might occur between p levels of the valence band and unoccupied s,d levels in the conduction band, provided of course 
that accessible final states lie above the vacuum level. They may give rise to intensity variations in the threshold 
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(a) 70 nm 



(b) 20 nm (c) 7 nm (d) 5 nm 

FIG. 10. Work Function maps for every thickness. 



-4.0 



U-3.9 



-3.8 o 



-3.7 



(e) 3.6 nm 




1-1.0. 



(a) 70 nm 



(b) 20 nm (c) 7 nm (d) 5 nm 

FIG. 11. MEM-LEEM transition maps for every thickness. 



- 1.6 - 



- 1.0 - 



(e) 3.6 nm 



spectra above cut-off energy and the shape of the rising edge of the photoemission threshold may be modified. In 
such a case the erfc parameters will no longer correctly describe threshold and inaccurate work function values may 
result. To check the effect of direct transitions on our work function values we took complementary image series 
using higher photon energy (Helium lamp hu = 21.2 eV and X-ray source A\-Ka hu = 1486.70 eV) for three of the 
BFO films: 20, 7 and 5 nm thin films, the thicknesses around the single to stripe domains transition. The higher the 
photon energy, the more only true secondary electrons contribute to the photoemission threshold. Results are similar 



within our energy resolutions (see Fig. 12 ) for 20 nm and 5 nm thin films. Notably, threshold widths for both types 
of sources are within 2%, showing weak influence of the direct transitions on the threshold spectra. In fact, it seems 
that p to s,d transitions in our BFO samples leave the measured position of the low energy cut-off in the spectra 
largely unchanged. Therefore, the influence of direct transitions on the work function can be neglected here. 



— I — I — I — I — I 1 1 



T 1 1 1 1 1 I I I 




0.4 



m 

-\ 0.2 
-\ 0.0 



FIG. 12. Thickness dependence of A^f (red squares for Hg lamp, green diamonds for X-ray source and Hel lamp) and AE 
(black circles). 



3D generalization of Landau-Ginzburg-Devonshire to BiFeOs Thin Films 

We start from the Ginzburg-Landau Free energy expressed in the form of the expansion with respect to the 
polarization P: 

F{P) = {Pi + P2) + ia,p2 + (P^ + P^) + i/^^P.^ + \p,PIpI + (P| + P^) P| + \lPt - {E, + E,) P, 

(3) 



If 



Py = then the equilibrium conditions result in the following equations: 



a,P. + li2Pl + 2/34P.P! + iPl = E, + Ed 



(4) 
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+ PiPl + /33P'P± + /34P±P' = (5) 

Where: 

Ed = ^^Q^^ ~ -^^^ _ Z (6) 

Here, h' is the width of the polarized region. The total thickness of the film is h = -\- where d is the dead layer 
width. We will assume that d <^ h. eo is the vacuum permittivity. Cd is the dielectric constant of the dead layer and 
65 is the so-called background dielectric constant (which is independent of the film thickness) [15] . U is the voltage 
between the contacts. Furthermore, Ed is the depolarizing field P^115j, and Ez = U/h. 

From ([5|, we have the choice, whether = or: 

Pi + P3 

This latter equality reveals that the z-component of the polarization influences the in-plane component, and vice 
versa, the magnitude of the in-plane component of the polarization influences the z-component. Now we substitute 
Eq. into Eq. Q and get: 



Where: 



Cdh + Cbd 



a = az ^ z — — ^ — 2p. 



d 

eo (erf/i' + Ebd) 
d 



2/3| 



a — OLz — 2p4 



Notice that ol^ is modified with respect to a^, and can even change sign, because of the depolarizing field and the 
correction due to the coupling of the in-plane component of polarization with its out-of plane component. Furthermore, 

is smaller compared to ^2 when all /3's are positive. This modification can even result in a negative and therefore 
change the second-order phase transition to a first-order one. In the case U = 0, Equation ([s]) has two stable solutions. 
One is = 0, while the other is: 



27 

One can easily show that such latter solution can be rewritten as: 



= AWB+Wl-^ (10) 



Where 
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D 

A 

B = 

^eff 



1 



(11) 
(12) 

(13) 
(14) 



Equation (10) is the one that has been used in the manuscript to fit the data of Fig. 5b. Note that in this fitting, we 
allowed A to take an arbitrary value, because Pmax in experiment is not very well defined (one cannot consider very 
thick films since they become too insulating for Photoemission Microscopy). However, the resulting A was numerically 
found to be very close to its ideal value provided in (12). Specifically, the ratio between the actual and ideal values 
for A was found to be about 1.05. 



Note that Equation (10) is, of course, valid provided that 



l-^>0 



B 



l-^>0 



These two conditions were met in the fit of the data of Fig. 5b for films thicker than 5.6 nm, since we numerically 
found that hefj = 5.6 nm and B = 0.16. 

It is also interesting to realize that the solution of Equation ([9| can adopt a more simple form than Equation (10) 
in some particular cases. For instance, if < 0, > and 7 = then: 



max 



Where 



9eff = T 

p2 ^ 

max 



Equation (15) has the same analytical form than the formula provided by Maksymovych et al [13]. However, the 



physical meaning of the parameters entering Equation (15) is different from those given in Ref. [13 , because, here. 



the polarization has three non-zero Cartesian components (rather than a single one). 



